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The thermal conductivity of a quasiparticle (QP) system described by a relativistic four-fermion 
interaction model in the presence of an external magnetic field is calculated. It is shown that, for 
narrow width of quasiparticles, the thermal conductivity, as a function of the applied magnetic field, 
exhibits a kink behavior at a critical field B c ~ T 2 . The kink is due to the opening of a gap in the 
QP spectrum at a critical magnetic field B c and to the enhancement of the transitions between the 
zeroth and first Landau levels. Possible applications of the results are discussed. 

PACS numbers: 74.25.Fy, 74.60.Ec, 74.72.-h, 11.30. Qc 



It is now well known from the study of field theoretical 
models that an external magnetic field can be a strong 
catalyst for a symmetry breaking leading to the gener- 
ation of a fermion dynamical mass even at the weakest 
attractive interactions jj| . This phenomenon of magnetic 
catalysis (MC) of symmetry breaking has been shown to 
be essentially model independent and to have wide appli- 
cations in several physical systems in particular, in 
condensed matter physics B-|^]. The essence of the MC 
effect lies in the dimensional reduction of the fermion 
pairing dynamics when the pairing energy is much less 
than the Landau gap y/eB [B is the magnitude of the 
magnetic field induction). In this case, the fermions are 
mostly confined to the lowest Landau level (LL) and their 
attractive interaction, regardless how small it may be, 
leads to mass generation. 

In recent years, several authors ||-||] suggested that 
the magnetic catalysis may play a relevant role in the 
physics of high-T c superconductors in external magnetic 
fields. This idea was inspired in part by the outcome 
of various experiments |9|~pd| initiated by Krishana's et 
al. on the thermal conductivity of high-T c compounds 
in the presence of an external magnetic field perpendic- 
ularly applied to the cuprate plane. According to these 
experiments, at temperatures significantly lower than the 
critical temperature of superconductivity, the profile of 
the thermal conductivity k versus the magnetic field dis- 
plays a sharp break in its slope at a transition field B c , 
followed then by a plateau region in which it ceases to 
change with increasing field up to the highest attainable 
fields - 14T. 

To understand how the magnetic catalysis was related 
[^-|| to these experimental observations, let us recall that 
high-T c superconducting cuprates are d-wave supercon- 
ductors which implies that their Fermi surface is charac- 
terized by nodal points where the gap function vanishes. 
Then, the low-energy quasiparticle (QP) dynamics is con- 
centrated mainly around these nodal points because these 
states can be occupied even at very low temperatures. 
The QP excitations in the vicinity of each node can be 
described by a massless Dirac Lagrangian [p~2|Jl3| . Hence, 



a natural candidate to modelling the QP interactions 
in a d-wave superconductor is a 3-dimensional relativis- 
tic four-fermion interaction (Nambu-Jona-Lasinio (NJL) 
type) model. It is well known [jl],[| that the gap equations 
of such a kind of models in the presence of an exter- 
nal magnetic field lead to the generation of a dynamical 
fermion mass (QP gap), which scales (at zero tempera- 
ture) with the field asm~ \feB. The critical tempera- 
ture at which the dynamical mass vanishes is determined 
by the dynamical mass at zero temperature, and there- 
fore, it scales with the magnetic field with the same law 
as the zero-temperature mass. One of the outcomes of 
the above mentioned experiments Q was that the critical 
temperature for the appearance of the kink-like behav- 
ior of the thermal conductivity scales with the magnetic 
field as T c ~ \fB, Then, the square root field depen- 
dence of the critical temperature in (2 + 1)— dimensional 
NJL models was considered as a hint that the MC phe- 
nomenon could be important to explain the experimental 
results of Krishana et al. 

By now, several possible mechanisms for the arising 
plateau have been proposed |j-g|,|l4] ^6| . All of them 
are based on the adoption of the QP picture as the low- 
energy excitations at the nodes of the d— wave symmetric 
order parameter. In the approach of Refs. follow- 
ing the ideas of Ref . [ p7|Jlq ] , it was assumed that the QP 
spectrum in the magnetic field was characterized by Lan- 
dau levels. Furthermore, in Ref. |i] ||,[l5| it was supposed 
that a QP gap opens at the critical field B c leading to 
the exponential vanishing of the quasiparticle contribu- 
tion in the thermal conductivity. Since the total con- 
ductivity is assumed to be the sum of the QP term kqp 
and the phonon term K p h with the phononic part inde- 
pendent of the field, this would lead to a plateau for 
high fields. Laughlin relates the QP gap to a weakly 
first order phase transition leading to the development of 
an additional small parity-violating superconducting or- 
der parameter. In another line of reasoning |p[-p|,^4| , the 
drop in the QP's conductivity, hence the plateau in the 
total conductivity, is ascribed to the opening of a field- 
induced energy gap in rf-wave superconductors due to a 
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second order phase transition. 

On the contrary, the approach of Franz fl(3| does not 
appeal to the generation of any QP gap. Albeit the 
phononic part of the thermal conductivity is also con- 
sidered to be field-independent, in Franz' scenario the 
QPs conductivity itself becomes field-independent above 
a crossover field B c . The effect is primarily due to the 
compensation of the enhancement of the QP density of 
states, arising from the Doppler shift of QP's spectrum 
in the superfluid velocity field around vortices (Volovik 
effect p!9[), and the growth of QP's width caused by the 
scattering on vortices. The Volovik effect certainly plays 
a decisive role at weak magnetic fields and low tempera- 
tures (less than 1 K), where the increase in the thermal 
conductivity was predicted theoretically [^o|, and later 
confirmed in the experiment p0| , to follow the \f~B de- 
pendence of the density of states. However, Volovik's 
arguments are based on semiclassical calculations valid 
in case one has well isolated vortices with supercurrents 
circulating around them. At higher fields, when vortices 
begin to overlap, this semiclassical picture should be re- 
placed by a more adequate quantum mechanical treat- 
ment. 

It should be stressed that none of the above works ad- 
dress the question of the kink itself in the thermal con- 
ductivity. Moreover, we should mention that the appli- 
cability of the Landau level quantization scheme used in 
the QP picture of the d-wave superconductor in a con- 
stant magnetic field has been recently subjected to in- 
tense criticism by several authors 
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The Landau 

level- like spectrum proposed by Gorkov and Schiffer [jl7| , 
and Anderson Jl8| was based on the assumption that the 
spatially dependent superfluid velocity can be neglected, 
but according to Melnikov the superfluidity veloc- 
ity strongly mixes individual Landau levels. In order to 
take into account the effects of spatially varying super- 
currents, Franz and Tesanovich ]22| ] introduced a singular 
gauge transformation to trade the phase of the order pa- 
rameter for a ficticious U(l) gauge field. As a result, 
the two fluxes, due to the external magnetic field and 
the U(l) field, cancel each other on average, washing out 
the whole picture of Landau level quantization. In this 
framework the low-energy quasiparticles are in fact Bloch 
waves described by massless Dirac fermions moving in a 
vector potential of physical supercurrents but zero aver- 
age magnetic field. While the arising scenario might seem 
promising, the approach of Ref. [p2p3| | itself is not free 
of problems: there is no good foundation from a physi- 
cal point of view to use singular gauge transformations 
and to introduce fictitious gauge fields, both of which are 
needed in this scheme. Moreover, because of the singu- 
lar transformations, the physical picture itself becomes 
strongly dependent on the way such transformations are 
implemented (i.e., different singular gauge transforma- 
tions can give rise to different spectra, depending on the 
rcgularization scheme J|j|). Finally, as far as we know, 
it has not been proved up to now that such a formula- 



tion can be applicable in the regime of weak-to-moderate 
magnetic fields, where the individual vortices overlap and 
the magnetic field can be considered to be fairly uniform. 

The aim of the present paper is to discuss a mechanism 
that generates the kink effect within the framework of the 
MC phenomenon. We emphasize that although, as shown 
in the derivations that follow, our numerical curves are in 
close qualitative agreement with the experimental curves 
of Ref. |J, we do not pretend to claim that our results 
explain the experimental findings of Krishana et. al. 
as our calculations are based on the Landau level quan- 
tization whose applicability to d-wave superconductivity, 
as mentioned above, is not well established. Our main re- 
sult in this paper is then of a more general and theoretical 
character. That is, we show that the MC phenomenon 
can be responsible for a kink effect in the thermal con- 
ductivity of (2+l)-dimensional relativistic fermion sys- 
tems, and we point out the very important finding that 
the kink-like effect close to the critical field is essentially 
model independent, since it is determined by the criti- 
cal behavior of the dynamically generated mass near the 
phase transition point. This fact makes our result rel- 
evant beyond the particular model under consideration, 
linking it to the universality class of theories with such a 
critical behavior. 

In connection with this, we would like to call the at- 
tention of the reader to a new class of recently discovered 
materials J24J], which could be a playground for applica- 
tions of the phenomenon of magnetic catalysis. The so- 
called highly oriented pyrolitic graphites (HOPG) have 
layered structure with two isolated points in the Brillouin 
zone where the dispersion is linear, so that the electronic 
states can be described here too in terms of relativis- 
tic Dirac particles [p5| . It was already suggested |2(| 
that the magnetic catalysis in the system of Coulomb- 
interacting Dirac fermions can provide an explanation 
for the semimetal-insulator phase transition observed in 
HOPG in the presence of a magnetic field perpendicu- 
lar to the layers. Since the quasiparticles in a graphite 
do not suffer from the Bloch waves versus Landau levels 
dilemma, our present results should have full strength 
there and we anticipate that the thermal conductivity of 
these systems should exhibit a behavior similar to the 
one we are reporting here. 

In this paper we present for the first time a consistent 
calculation of the thermal conductivity in the presence 
of an external magnetic field in a model with the sim- 
plest four-fermion interaction for quasiparticles. We use 
the same constant magnetic field approximation that was 
already explored in Ref. |(|[?J to calculate the thermal 
conductivity. However, our calculation deviates consid- 
erably from what was done in ^,0. Not only we take 
into account the contribution of all Landau levels, but 
the definition of the heat current itself is different. In 
the limit of narrow width of quasiparticles (r <C T, \f~B), 
after a gap is opened, the thermal conductivity exhibits 
a new term proportional to a 2 {a is the gap). This term 
originates from the compensation of the matrix elements 
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of transitions between the zeroth and the first LLs (be- 
having as 1/B) and the LL density of states which in 
turn is proportional to B. This is one more manifesta- 
tion of the MC phenomenon: not only a gap is induced 
by the magnetic field, but the transitions between the 
zeroth and the first LLs are enhanced. In mean-field ap- 
proximation and near the phase transition point the gap 
behaves like a ~ 0.523\/eB — eB c what yields a positive 
contribution into the slope of the thermal conductivity 
leading to a jump (kink effect) of n(B) at B = B c . 

We start from the following NJL model in an external 
magnetic field in (2 + 1) dimensions, 
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In (|l|) the vector potential is taken in the symmetric 
gauge = (0, — -jXi) , and is a characteristic 
velocity P7| . In what follows we take K = Vd = ks = 1 
and absorb c into the charge e. We will restore these con- 
stants when needed. We assume also that the fermions 
carry an additional flavor index a — 1, . . . ,N (N = 2 for 
realistic d-wave superconductors). The Dirac 7 matrices 
are taken in a reducible four-component representation. 

The Lagrangian density ([!]) is invariant under the dis- 
crete (chiral) symmetry tp — > 75^, tp — > — "075 7 which for- 
bids the fermion mass generation in perturbation theory. 
The mass generation can be studied introducing an aux- 
iliary field a — —(g/Ntyiipi by means of the Hubbard- 
Stratanovich trick which permits one to integrate over 
fermion fields in the path integral representation of the 
partition function. The field a has no dynamics at the 
tree level but it acquires the kinetic term due to fermion 
loops. Likewise, a nontrivial minimum of the effective po- 
tential (the expectation value of a) gives mass to fermions 
and spontaneously breaks the discrete symmetry lead- 
ing to a neutral condensate of fermion- antifermion pairs. 
Studying the minimum of the effective potential we find 
that, at a fixed temperature T, there is a critical value of 
the magnetic field \/eB c /T ~ 4.148 such that for subcrit- 
ical fields eB < eB c the gap is zero, while for eB > eB c 
a nontrivial gap appears. The critical curve equation has 
the form (v D /c) 2 10 1Q B = 21.5 ■ T 2 for B measured in 
Tesla. Using the approximated value of the characteris- 
tic velocity vd — 10 7 cm/ s |§,^6), we obtain the critical 
curve B = 0.014 • T 2 . Notice that the obtained critical 
curve fits well the experimental curve of Ref. || . 

To derive an expression for the static thermal conduc- 
tivity in an isotropic system, we follow the familiar linear 
response method and apply Kubo's formula P8| 
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where V is the volume of the system and Ui(x,t) is the 
heat-current density operator. The brackets denote aver- 
aging in the canonical ensemble. Physically, the thermal 



conductivity k appears as a coefficient in the equation re- 
lating the heat current to the temperature gradient u = 
— kWT under the condition of absence of particle flow. If 
we neglect the chemical potential, the heat density coin- 
cides with the energy density e, hence the quantity u that 
satisfies the continuity equation e(x) + V .u(x) = can 
be interpreted as the heat current. From the Lagrangian 
density (|l|) we find Ui = | (^idpip — do^pji^p) . 

Neglecting vertex corrections p9| the calculation of 
the thermal conductivity reduces to the evaluation of the 
bubble diagram ]30| ] . In this approximation, making use 
of the spectral representation for the fermion Green's 
function, we arrive at the following expression for the 
thermal conductivity in the Matsubara formalism 
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The spectral function A(w, k) = — 7r _1 ImS' i? (a; + is, k) is 
derived from the fermion propagator in a magnetic field 
decomposed over Landau level poles Jl|,^l|, 
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n— 



{^M n -a)h{k)-f 2 (k) 



(r/°Mn + (T)fl(k)+f 2 (k) 

(w - M„) 2 + F 2 
, M n = Vcr 2 + 2eBn, (4) 



(cu + M n ) 2 +T 2 
h(k) = 2 [P-L n (s) - P + L n -!(a)] ,f 2 (k) = 4(fc-7)L I \_ 1 ( s ) 

Here P± = (1 ± ij 1 j 2 )/2 are projector operators, L n , L\ 
are Laguerre's polynomials (L 1 _ 1 = 0), s = 2k 2 /eB, and 
a is the dynamical fermion mass obtained from the finite 
temperature gap equation in the magnetic field. We in- 
troduced also the quasiparticles width replacing e in the 
definition of A(u>, k) by finite V, which is due to inter- 
action processes, in particular, scatterings on impurities. 
Straightforward calculation of the trace and further in- 
tegration over momenta in Eq. ([}]) produce Kronecker's 
delta symbols 5 n , m -i+5 OT ,n-i due to orthogonality of La- 
guerre's polynomials. This implies that only those tran- 
sitions between neighbor Landau levels contribute to the 
heat transfer. After performing one of the sums in (|^) 
we obtain 
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Further summation over n in Eq. (^) can be carried 
out expanding the integrand in terms of partial fractions. 
The resulting sums are expressed through digamma func- 
tions and the final expression for k is written as follows 
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This formula is the main result of our paper and we can 
use it now to study the different asymptotic regimes. It 
is easy to show that in the limit of zero field k reduces to 
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where we put a = 0, since the mass is not generated 
in the zero-field weak coupling NJL model. Eq. (|t]), 
up to an overall factor k 2 B {v 2 F + v^/Hvfvai coincides 
with the corresponding expression obtained in Ref. [ p^[ . 
With the overall factor replacing N (= 2 in real d— wave 
superconductor), Eq. (|t]) reproduces the universal (or 
residual) thermal conductivity at low T in the so-called 
"dirty" limit, T < Y @. The residual conductivity 
was recently observed in experiments [ |32"| giving explicit 
confirmation of the existence of gapless quasiparticles in 
d-wetve cuprates at T < T c . 

Let us consider the case F <C T with B ^ 0. In this 
approximation we obtain 
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Asymptotically, at V eB 3> y 'eB c ~ T, the dynamical 
mass behaves as a ~ \fe~B what leads to an exponential 
decrease as in the case of gapless fermions. However, 
the term cosh~ 2 er/2T in Eq. (||) is of order one when 
a < 2T, thus, there is no suppression of this term for a 
certain range of fields where it is almost field independent 
(plateau region). 

Near the phase transition point a ~ a\/eB — eB c and 
the first term in Eq. (|J) gives positive contribution to the 
slope of the thermal conductivity at eB > eB c leading to 
the jump in the slope of n (kink effect) at eB = eB c . The 
parameter a is model-dependent, but the scaling of a is 
determined by the universality class of the theory. For 
the NJL model (0) we find, in mean-field approximation, 
a - 0.523. 

The explicit appearance of the kink has been corrobo- 
rated by numerical calculations as shown in Fig.|l|. Notice 
the break in the slope of k (kink effect) at the critical 
value B c in the presence of a. For B > B c the kink is 
followed by a region where k is only weakly dependent 
on the field. While decreasing the temperature, the posi- 
tion of the kink moves to the left in accordance with the 
critical line B c — 0.014T 2 . Notice the similarity of our 
results (Fig. Q) with the experimental behavior reported 



in H for high-Tc superconductors. Moreover, Fig. [l] 
also reproduces the experimentally observed crossing of 
the curves that occurs with decreasing T in such a way 
that the lower T curve reaches the higher value at large 
fields. 

Because of the important role played by the first term 
in (||) to get the kink effect, it is instructive to clarify its 
origin. From Eq. (§) one can see that while the contri- 
bution of transitions between Landau levels with n > 1 
in the integrand behaves as l/(eB) 2 for large fields, the 
contribution of transitions between the zeroth and the 
first LLs decreases only as the first power of the field 
(~ 1/eB). Since the density of LLs is proportional to 
eB, this implies that the transitions between the zeroth 
and the first LLs are not suppressed, even though the 
gap between the levels grows with the field. Note that 
for gapless QPs the transitions with n = are completely 
suppressed in the regime Y — > 0: indeed, their contri- 
bution is proportional to 8(u>) and the integrand in (|^) 
contains an uj 2 factor. For gapped QPs, S(lu) is replaced 
by 5{uj ± a) and the transitions with n = survive. 

In summary, in the present paper the thermal con- 
ductivity at B ^ in the framework of a relativistic 
four-fermion model was calculated. Assuming a uniform 
magnetic field approximation, we showed that the ther- 
mal conductivity exhibits a kink behavior when the field 
reaches the critical value B c = 0.014 • T 2 . Two main 
features determine this effect: the generation of a QP 
gap in a magnetic field (MC phenomenon) and the lack 
of suppression of zeroth-first LLs transitions. We stress 
that the appearance of a kink in the thermal conductiv- 
ity is model independent, being only determined by the 
critical behavior of the induced dynamical mass near the 
phase transition. The universality character of this result 
leaves an open window for possible future applications. 
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FIG. 1. The magnetic field dependence of k at T — 20K 
and T — 15K in the narrow width case (r = 5K). The 
solid lines represent k/ko when a QP gap a is MC-induced at 
B > B C (T) (B c (20) = 5.75T, B c (15) = 3.23T). The dashed 
lines represent the behavior of k/ko when a remains zero at 
B > B C (T). 
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